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We re-examine the tunneling data on A15 superconductors by performing a generalized McMillan- 
Rowell tunneling inversion that incorporates a nonconstant electronic density of states obtained 
from band-structure calculations. For NbsSn, we find that the fit to the experimental data can 
be slightly improved by taking into account the sharp structure in the density of states, but it is 
likely that such an analysis alone is not enough to completely explain the superconducting tunneling 
characteristics of this material. Nevertheless, the extracted Eliashberg function displays a number 
of features expected to be present for the highest quality NbsSn samples. 
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I. INTRODUCTION 

Twenty years ago, the A15 superconductors in the A3B 
structure, with A a transition metal and B an sp-metal, 
were the highest transition temperature superconductors 
known. The first such superconductor found was VaSi, 
which was discoverecB in the 1950's to have a Tc of about 
17 K. In the ensuing years, other equilibrium compounds 
like VsGa and NbaSn were discovered to have Tc's in the 
range of 15-18 K. B-element-poor compounds like NbsAl, 
NbsGa, NbsGe, and NbaSi were also found, with Tc's as 
high as 23 K. The structural, electronic, magnetic, elas- 
tic, vibrational, and superconducting properties of these 
materials were widely studied .ua 

The use of tunneling spectroscopy to probe supercon- 
ducting properties of the A15 materials was hindered by 
a number of materials issues. The fabrication of high- 
quality tunnel junctions was difficult because the use 
of native oxides for the tunneling barrier did not yield 
reproducible results, and hence artificial barrier layeis 
had to be grown on top of the A15 superconductorsoEl 
As the quality of the tunneling data improved, it be- 
came clear that these materials do not exhibit the sim- 
ple behavior seen in conventional strong-coupling s-wave 
superconductors like Pb, Hg, and Nb. In particular, 
the reduced tunneling density of states displays a rapid 
overswingQ followed by a sharp return to zero at ener- 
gies near and above the maximum phonon energy- A 
revision of the McMillan-Rowell tunneling analysisQ by 
Arnold and Wolfl allowed this data to be fit to high pre- 
cision by assuming the presence of an additional normal- 
metal layer (characterized by a width and a scattering 
strength) between the superconductor and the insulat- 
ing barder. However, to fit the most recent experimen- 
tal datad within this proximity-effect modified tunneling 
theory, one needs to assume that the width of the nor- 
mal region approaches zero while its scattering strength 
becomes unusually large. 



Alternatively, it has been suggested that the rapid 
overswing in the high-energy regime in the tunneling den- 
sity of states is related to the presence of sharp struc- 
tures-in the electronic density of states near the Fermi 
levelfl Band-structure calculationsEI show that the elec- 
tronic density of states in these materials has peaks of 
width on the order of 100 meV near the Fermi energy. 
Such sharp structures require, at tlie very least, a refor- 
mulation of the Migdal-EliashbergiJ many-body analysis 
to incjkjd&i effects of a nonconstant electronic density of 
statesll2l"ll3 within an energy range on the order of the 
maximal phonon energies of the material. 

Other theories of superconductivity in A15 compounds 
go beyond simply generalizing the standard Migdal- 
Eliashberg theory to allow for energy dependence p. the 
density of states. For example, Yu and AndersonEil ex- 
amine what happens in an electron-phonon system that 
is coupled strongly enough to have the single-electron 
(polaronic) phase become unstable to bipolaronic (pre- 
formed pair) phases. Such systems display quite different 
behavior, but these theories have not been developed to 
the point where direct comparison to materials-specific 
tunneling conductances can be made. 

Since the discovery of the high- Tc cuprates, work on 
the A15 compounds has virtually ceased - this despite 
the many fundamental questions about these materials 
that remain open. In this paper, we investigate whether 
high-quality tunneling data for the A15 materials can 
be understood within a conventional Migdal-Eliashberg 
framework generalized to include an energy-dependent 
electronic density of states. We perform this analysis 
and extract an experimentally fit Eliashberg function 
Q!^T(0), and Coulomb pseudopotential /x*, for NbaSn, 
which has Tc = 18 K. Data from high-quality tunnel 
junctions grown on this material are availablefl 

In Section II, we derive the formalism needed to numer- 
ically perform the tunneling inversion including a non- 
constant electronic density of states. This derivation and 
the computational algorithm are different from those de- 
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veloped previously in that an exact analytic continuation 
that properly treats the Coulomb pseudopotential and 
allows the calculations to be performed relative to the 
normal state is employed. In Section III, we present our 
numerical results for the tunneling inversion with both 
constant and nonconstant electronic density of states. 
Conclusions are presented in Section IV. 



II. FORMALISM 

Our analysis begins with the calculation of the band 
structure and electronic density of states (DOS). jCal- 
culations are carried out using the VASP package|l3 a 
plane-wave basecL-density-functional code using ultrasoft 
pseudopotentialsJij The electronic wave functions are ex- 
panded in plane waves up to a cutoff of 219 eV, and 
the electron-electron interaction is treated within the 
local density approximation (LDA) witi the Ceperley- 
Alder exchange-correlation functionalO The Brillouin 
zone is sampled on a Monkhorst-Pack meshlla of at least 
20 X 20 X 20 points. We find the optimized crystal struc- 
ture to bcptetragonal, with a small distortion of the Nb 
sublattice.Ej The peaks in the density of states near the 
Fermi level are very sensitive to this sublattice-distortion. 

In conventional Migdal-Eliashberg theoryJi3 the elec- 
tronic density of states is chosen to be a constant (with 
an infinite "bandwidth"), and the energy cutoff of the 
theory is provided by the finite range of the Eliashberg 
function a^F{ fl), which measures the ability of a phonon 
of energy to scatter electrons on the Fermi surface. 
When performing the many-body theory calculations, 
one begins on the imaginary axis,|-where the Coulomb 
pseudopotential has a sharp cutoff,c2l andj-then performs 
an exact numerical analytic continuationEd to calculate 
real-axis properties. This technique allows for a proper 
treatment of the "soft" cutoff for the Coulomb pseudopo- 
tential on the real axis. In addition, because the su- 
perconducting and normal states appear more and more 
similar at high energies, and because of the finite fre- 
quency cutoff (chosen to be 6f2max here), the self-energy 
begins to deviate from the exact result in the normal 
state as the energy increases. One can correct for this 
by performing the perturbative expansion relative to the 
normal state. In this case, one adds the exact normal- 
state self-energy minus the normal-state self-energy cal- 
culated with-the energy cutoff used in the superconduct- 
ing theory.cS Such a scheme-jwas used when examining 



effects of vertex correctionsc^ in Pb. When an energy- 
dependent electronic density of states is used, such a com- 
putational scheme becomes problematic because there is 
no longer an exact analytic expression for the normal- 
state self-energy. Furthermore, a choice must be made 
for the energy cutoff, because, assuming a bandwidth on 
the order of a few eV and a temperature on the order of 
0.1 meV, the number of Matsubara frequencies required 
for an energy cutoff set by the bandwidth would be too 



large to perform calculations efficiently. 

We adopt an alternative scheme here. We start by 
first calculating the normal-state self-energy at T = 0. 
Since the Matsubara frequencies become a continuum at 
T — 0, such a calculation can be performed simply by re- 
placing the Matsubara summations by integrals along the 
imaginary axis, which are computed using conventional 
quadrature techniques. Next, we choose the energy cut- 
off to be the same as that used in the constant density of 
states calculations, namely dflmax- We calculate the self- 
energy in the normal state at T (relative to the normal 
state at T = 0) with the finite energy cutoff, and then 
add the T — normal-state self-energy (with no cutoff) 
to the finite-temperature self-energy. Then we calculate 
the superconducting self-energy relative to the normal 
state at T using the same cutoff. Finally, we add the full 
normal-state self-energy at T to get the superconducting 
self-energy. 

Our many-body formalism on the imaginary axis fol- 
lows most closely to Ref. |l3[ We use a Nambu-Gor'kov 
formalism and evaluate the self-consistent perturbation 
theory (using dressed phonons) in the Hartree-Fock ap- 
proximation. (The Hartree term provides just a constant 
shift to the chemical potential and is ignored.) Pertur- 
bative calculations are performed for the normal state 
and the superconducting state. We begin our analysis 
in the normal state, where the resulting self-consistent 
equations are 



Jde 



n 



fJ-N 



P{Ef) uj^Z^ + 



X'nf 



(1) 



zL 



Jde 



T 



^nZ'^ 



p{Ep) u:lZ'^ + {p^ 



(2) 



A prime indicates the normal-state perturbation theory, 
iijjn — i7TT{2n + 1) is the Fermionic Matsubara fre- 
quency, x'm — x'(*'^m) = ReE'(iix>m) is the real part 
of the normal-state self-energy, and Z'^ = Z'{iujm) = 
1 — lmT,'{iujm)/'-^m is the so-called renormalization func- 
tion determined from the imaginary part of the normal- 
state self-energy. The symbols Xm-n are the electron- 
phonon coupling strengths 



Am-n — X{iuj,n — ioJn) 



dfla^Fin) — 

+ [LO: 



2n 



(3) 



and A = Aq. The function p(e) is the electronic density of 
states as determined by the band-structure calculation, 
and p{Ep) is the density of states at the T — Q chemical 
potential of the band-structure calculation [Ep). The 
chemical potential in the normal state is /ijv. In the 
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limit as T — ^- 0, the number of Matsubara frequencies 
becomes an infinite continuum, and the summations can 
be replaced by integrals. We find 



'I- \ 1 



dui' X(iLU 



IJ-NO 



p{Ef) uj'^Z'^ii^') + [iiNO - e - x'oii^'W 



(4) 



ZTTU) J 



Jde 



-ILV') 



p{Ef) u'^Z'^iiLo') + [/iATo - e - x'S^'W 



(5) 



with the subscript zero denoting the results are at T = 
0. Evaluating these integrals with quadrature routines is 
much more efficient than calculating the Matsubara sums 
with a large energy cutoff at finite temperature. 

Our strategy for determining the functions Xq {iu) and 
2^(10!) is to create a nonuniformly spaced grid on the 
imaginary axis (we use 336 points) , with an upper cutoff 
many times the electronic bandwidth. The grid is con- 
structed in the following fashion. The first grid point 
is chosen to lie at iu)o = 0. Subsequent grid points are 
chosen by adding the new step size to the old grid point 
iu>j+i = iijJj + idj, with 6j increasing by a factor of 1.1 
with each step from its initial value So = 0.09 meV [hence 



Sj = 0.09x(l.lp]. Werestrict^j < 400 meV- 



-if Sj would 



be larger than 400 meV, we set it equal to 400 meV. The 
grid is reflected about iui = to construct the negative 
frequency axis. The integrals are then evaluated using 
a simple Riemann sum over the nonuniform grid. Once 
the functions Xo ■^o known on the grid points, 
we linearly interpolate to evaluate them at any point on 
the imaginary axis. 

We employ the T = normal-state solutions as an ap- 
proximate solution for high frequency. This is included 
by first calculating the normal-state perturbation the- 
ory at finite temperature, using the same Matsubara 
frequency cutoff as used in the superconducting state 



LOnZQ^iUn) 

'^lZ'^{iuJn) + {PNO - e - 
+ [Z^iitOm) - 1]. 



(7) 



These equations would be exact if the summations over 
Matsubara frequencies for the T = quantities were re- 
placed by integrals. Since the normal-state self-energy 
does not depend too strongly on T for low temperature, 
this approximation is accurate for low T. Note that the 
chemical potential /xjv typically changes by about one 
meV (at T = 1.2 K) from the zero-temperature value 
/Ujvo (in the nonconstant density of states case). 

The final set of equations we need are for the super- 
conducting phase. The self-consistent equations are cal- 
culated "relative to the normal state" at T: 



An 



de 



MS 



Xn 



^IZl (/i5 - e - Xn) 

PN-e-x'n 
^lZ'^ + [p^-e-x'n? 



"I" Xmi 



(8) 



T - f 

, ^ J 



de 



P{Ef) 



^nZn 



u;lZl + {ps-e-Xn? + KZl 

^nZ' 



+ [mat - e - x!nY 



+ \K - 1], 



(9) 



and 



in as usea m me superconaucimg siaie /■ 
and then adding the difference between ^mZm = T^ (Am-n— /^*) / 



the r = normal state solution and the finite-T solu- 
tion to the superconducting solution, as shown below. 

We solve the following self- consistent equations for the 
normal-state self-energy at temperature T: 



\u>n\<Mlc 



p(Ef) 



^nZn 



^lzl + {^is-e-xn? + i^lzl■ 



(10) 



T 



E 

\uJn I <Wc 



A. 



de 



p{Ef) 



Pn 



Xn 



colZ^^ + ip^-e- 
Pno — e ■ 



■X'n? 
Xoii^n) 



UjlZ'^{iUn) + [mjvo - e - Xoii'^nW 



(6) 



zL = i + 



E 

|w„|<Wc 



de 



p{Ef) 



Here = A(ia;m) = Y,\2{ii^m)/Zm is the supercon- 
ducting gap determined from the off-diagonal self-energy, 
p* is the Coulomb pseudopotential, and ps is the chemi- 
cal potential in the superconducting state. Note that we 
add and subtract the normal-state results at finite tem- 
perature in order to ensure that the Matsubara frequency 
summations converge rapidly (wc = Sflmax is the cutoff 
frequency; the high-frequency tails are already included 
in the normal-state self-energy). This procedure allows 
for a rapid computation of the many-body Green's func- 
tions when there is an energy-dependent electronic den- 
sity of states. We only calculate the Green's functions at 
the Matsubara frequencies here. 
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The chemical potentials in the normal and supercon- 
ducting states are determined by the requirement that 
the electron density be equal to the equilibrium density 
of electrons for the given band. Our first step is to find 
the T = Fermi level for the band-structure density of 
states, which satisfies 



and 



(15) 



(16) 



dep{e), 



(11) 



with the factor of 2 coming from spin. In the normal 
state at T = 0, we use the T ^ limit of the identity 
Pe = 2T^^G(iti;„) and Eq. (|ll|) to produce the self- 
consistent equation for pNo- 



Ep 

2 / dep{e) 

' P-NO 



1 



dep{e) I duj 



1 



1 



iuj + UNO — e 



(12) 



The normal-state chemical potential at finite T (we use 
T = 1.2 K for the tunneling inversion) is found by 
comparison with the normal state at T = 0: = 
2rX;„[G(iw„) - G{iuJn)\T=o] (note the T = sum is 
an approximation to the continuum integral). Finally, 
for the superconducting state, we use the comparison of 
the normal-state filling with the superconducting-state 
filling to find p.s- 



\uJm I <LJc 

M5 - e - Xn 



Mat - £ + Xrn 
^liZ'^ + {PN ~e + X'm}" 



(13) 



The next step is to calculate the self-energy on thfenreal 
axis using an exact analytic continuation techniques We 
begin with the normal state at T = 0. The self-energy 
satisfies 



de 



2ttp{Ei. 



duj' \(uj — iuj') 

oo 

Pie) 



1 



p{Ef 



iLu'ZQ{iuj') + pNO 



dQa^F{Q) / de 



X'o( 



Pie) 



where the * denotes complex conjugation. The term 
A(cl' — iLu') is found from the spectral formula and the 
given Eliashberg function a^F{n): 



X{u! — iuj') 



dVL- 



a^F{n) 



(17) 



Now we calculate the finite-T normal-state results rel- 
ative to the T = calculation: 



X de 



T 



PiEf 



A(lJ - iuJn) 



Pie) 



iUJnZ^ + PN - e - Xn 

Pie) ' 

iw„Z^(iw„) -I- pNa - e- Xoi^^^n) 

1 

dna^F{n) 



MEf) 



tanh 



LJ — 

2T 



coth ( — I 
1 



depie) 



{uj — fl)Z'{LU — Q) + pn — e — x'i^ — il.) + irj 
1 

{lu - n)Zl^{iLj - ft) + pno - e — x'oi^ - ft) + iri 

Ki^ + ^v), 



where the Z and x functions on the real axis are deter- 
mined from equations analogous to Eqs. ([l5| ) and (^6|). 

For the superconducting state, calculations are per- 
formed relative to the normal state at T: 

II {u! + irj) — T ^ X{uj — iujn) 

\^^n I <Wc 

iuj„Zn + PS ~ Xn 



(18) 



X / dep{e) 



ulZl + {ps-e- x„)2 -f A2 Z2 

iuJnZ'^ + PN -X'n 



'^n^n +iPN -e- x!n? 



+ 



1 



MEf) 



dna^F{n) 



{lu - i})ZQ{uj — fl) + PNO - e - x'ai^ — ^) + iv' 

where 77 — > O"*". This is a self-consistent equation, because 
the second integral contains the self-energy on the real 
axis from the definitions 



(14) 



tanh 



n 



coth 



- n 



depie) 



2T J \ 2T 

(uj — il)Z{LU — ^l) + PS — e — xi^ ^ ^) 

(w - n)^z^{uj -n) + [ps -e-xito- fi)]^ + A{uj - n)^z{ijj - n 
(w - n)z'(ui - n) + PN - e - x'(w - n) 

(uj - f])2Z'2(u; -n) + [pN-e~ x'i^ - ^)? 
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+ T,'{Lj + ir]), 

for the diagonal self-energy and 
'Ei2{uj + if]) ~ T ^ [X{uj 

X / depie) 



1 



2/3(£;f) 
tanh 



^IZ^ + (MS - e - Xn)' + A2 Z2 

cj — r2 



coth 



2r 



rfep(e) 



2r 

X A(w - r2)Z(w - f7) 

for the off-diagonal self-energy. 

The tunneling conductance satisfies 



(20) 



[dl/dV], 
[dl/dV 



N 



■(uj) = Re 



^t^2_ A2(t^) 



(21) 



The reduced density of states (RDOS) is the ratio of the 
tunneling conductance to the BCS tunneling conductance 
minus one, which becomes 



RDOS{lo) = Re 



^0 



^Cj2 - A2(tj) 



(22) 



where Aq is the superconducting gap, defined as the so- 
lution to the equation ReA(w) = uj {a.t lo — Aq). 

To solve all of these self-consistent equations and to 
invert the tunneling conductance to get the Eliashberg 
function, we begin with the LDA band-structure elec- 
tronic density of states as input. We take the band to 
consist of the states that contain nine electrons below 
the Fermi level and four electrons above, corresponding 
to approximately ±0.8 eV about Ep. Hence the electron 
filling satisfies pe = 9. Our first step (i) is to determine 
the band-structure Fermi level Ep at T = by solv- 
ing Eq. (|ll|). Then we take two initial guesses for ^at. 
(ii) For each ^jv, we self-consistently solve for Eqs. (|l|) 
and (H) at T = 0. Once those functions are known, (iii ) 
we use a one-dimensional root-finder to solve Eq. ( p^ ) 
to determine ^n- The next step is to determine the su- 
perconducting chemical potential jis s-nd the Coulomb 
pseudopotential /i* . We first choose a guess for p.* . With 
/i* fixed, we choose two initial guesses for /is and (iv) 
self-consistently solve Eqs. (||), (^), and ([l^). Once those 
functions are determined, (v) we use a one-dimensional 
root-finder to solve Eq. ( [l3| ) and find /is. (vi) Now the 
real-axis equations for both the normal state at T = 0, 
Eqs. ( p^ ) and (|l^), and the normal and superconduct- 
ing states at T, Eqs. (18), (|l^), and (po|), are solved self- 
consistently. Once the gap function is known, we 
can find the superconducting gap Ao. (vii) Steps (iv), 
(v), and (vi) are repeated for different values of /i* un- 
til we find the particular value where Aq — 3.15 meV, 



which is th^MI)perconducting gap for NbsSn at low tem- 
perature. If we knew the correct Eliashberg function, 
then steps (i)-(vii) would be all that we need to calcu- 
late the RDOS from Eq. dH). But we need to find the 
best fit a^F{^) that is consistent with the experimen- 
tal RDOS. To do this, we simply employ the McMillan- 
Rowell procedure for the tunneling inversion. We start 
with a guess for a^F. (viii) Next we calculate the func- 
tional derivative of how changes in the Eliashberg func- 
tion affect the RDOS. The functional derivative is found 
by adding a small-weight Gaussian to a^F{^) centered at 
a given frequency Oq, repeating steps (ii-vii) to determine 
the RDOS, and calculating the functional derivative by 
taking the difference of the new RDOS with the RDOS 
for the original Eliashberg function. This step (viii) is 
repeated for each frequency in the discrete a^F to deter- 
mine the functional derivative matrix. The next step (ix) 
is to choose the weights for adjusting a^F to fit the ex- 
perimental RDOS better. Since the functional derivative 
matrix may be rank-deficient, we use a singular- value de- 
composition to determine the updated weights. Once the 
updated a^F{Q.) is found, (x) we apply an exploratory- 
data-analysis robust smoother followed by a Hanning 
smoother. Steps (viii-x) are repeated until the updated 
a2_F(Sl) ceases to produce a RDOS that is closer to ex- 
periment. We force the Eliashberg function to be positive 
everywhere and to be quadratic for < 2.5 mcV. The 
step size for the discrete frequencies at which a^F is eval- 
uated is 0.45 meV. The algorithm is depicted pictorially 
in Fig. I 

There is a numerical problem with this algorithm. In 
general, the weights for the shift in a^F are strongly os- 
cillatory, which forces the Eliashberg function to have 
large-amplitude narrow-width peaks. Since we do not 
expect the Eliashberg function to have such sharp peaks, 
they need to be smoothed away in the next update. Be- 
cause the functional derivatives are not calculated with 
this smoothing procedure (and there is no obvious way to 
incorporate the smoothing into the derivatives), we are 
limited in how closely we can reproduce the experimental 
RDOS, even though we have the full freedom to adjust 
a'F. 
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Eq. (9) 



(a) 



Eqs. (1) and (2) 



(vii) 



(ii) 



(iii) 



(iv) 



XQ(i(i))ZJj(i(o) 



Eq.(10) 



Steps (iv) 




(v) and (vi) 












X(ffl) Z(ffl) A((o) 



(Vi) 



Eqs. (6-i 



Eqs. (13) (14) 
(16) and (17) 



Eq. (11) 



(V) 



(b) 



(vii-b) 







a F(n)+ g(a-no) - 




(ii-vii) 


A((o) n* RDOS(<o) 




compute 




functional 




derivative 



a F (n) 
new 



(X) 



smooth 



(vii-c) 



S RDOS (to) 
5a^F(n) 



Singular Value Decomposition 



8a^F(n) 



(ix) 



FIG. 1. Block diagrams of the algorithm employed in the 
modified McMillan- Rowell tunneling inversion with a noncon- 
stant electronic density of states. In panel (a), we show the 
main algorithm used to determine the RDOS. Steps (i-vii) are 
described in detail in the text. Panel (b) shows the McMil- 
lan-Rowell strategy for updating the Eliashberg function in 
the tunneling inversion, with details described in the text. 



III. RESULTS 

We first consider a test case, where the experimental 
tunneling data for lead is used to perform both a conven- 
tional and a nonconstant DOS tunnehng inversion. In 
both cases we use the same code to perform the analysis 
as summarized in Fig. 0. In the constant DOS case, we 
use the DOS for lead (2.5 x 10^'* states per spin per meV 
per unit cell) and choose a band that contains two total 
electrons (including spin) and pe = 1. In the nonconstant 
DOS case, we use the electronic DOS for NbaSn (see in- 
set to Fig. I). This is of course an artificial problem that 
is presented for illustrative purposes. 

The results for the best fit Eliashberg functions are 
given in Fig. ^. The solid line results from a tunneling 
inversion done with a constant DOS, and is consistent 
with earlier analyses of the lead tunneling data (A — 1.56, 
maximal error of 0.002, and a root-mcan-square error 
of 0.0004). The tunneling-inversion result based on the 



energy-dependent NbaSn DOS is plotted as a dotted line, 
and corresponds to A = 1.42, with a maximal error of 
0.004, and a root-mean-square error of 0.0015. The fit 
with a nonconstant DOS is about four times worse than 
with the constant DOS. Note further, that the main effect 
of the sharp peak in the DOS is to reduce the overall 
scale of a^F for nearly all frequencies except the highest, 
where it is strongly enhanced. Qualitatively, the two a^F 
curves are very similar in shape. 




Q [meV] 



Fig. 2, Freericks et al., Phys. Rev. B 

FIG. 2. Eliashberg functions for a tunneling inversion in 
lead with (solid line) constant density of states and (dotted 
lines) the nonconstant density of states of NbsSn. The general 
shapes of the two curves are similar. The values of A are 1.56 
and 1.42 for the constant and nonconstant density of states 
cases, respectively. 

The quality of the tunneling data for lead is so good 
that one can actually see that the tunneling inversion 
works better with a constant DOS than with a noncon- 
stant DOS. What is interesting is that one can get rea- 
sonable results using the wrong nonconstant DOS due to 
the freedom allowed in choosing the function a^F. 

We next turn to the results for the tunneling inver- 
sion of NbsSn. The extracted Eliashberg functions are 
shown in Fig. |^. The thick curve shows the noncon- 
stant DOS results, while the thin curve is for the con- 
stant DOS case. The properties of these electron-phonon 
spectral functions are summarized in Table I. In contrast 
to the results for lead, where the Eliashberg function is 
reduced in the nonconstant DOS analysis, it is enhanced 
for NbaSn. The main difference, aside from an overall 
scale factor, is a large enhancement for > 20 meV, 
particularly in the highest energies, where the constant 
DOS a'^F is strongly suppressed. Note that the maximal 
allowed phonon frequency is the DFT-calculated 38 meV 
in both cases; the fitting procedure sharply suppresses 
the a^F at high energy for the constant DOS calculation. 
The general shapes of these curves are similar to those 
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found previouslyO but the overall scale is larger. The en- 
hancement occurs because we are unable to reproduce 
the full overswing. In order to compensate for this, A is 
increased and the quality of the low-energy fit is reduced. 
Note that we fit all experimental points, including those 
in the overswing region {to > 35 meV) — nevertheless, we 
arc unable to get "good" results for the extracted a^F 
in the sense that the best fit a^F is unable to accurately 
reproduce the tunneling data over the entire experimen- 
tally measured range (due to the inability to produce the 
full overswing). 

There are, nevertheless, a number of promising fea- 
tures of this calculation. If we compare Fig. ^ to the 
data for low- Tc (disordered) NbsSn junctions, we find the 
broad low-energy peak and shoulder (present from 4 to 
8 meV) lies a few meV lower in the high- Tc material than 
the low-Tc oneB (where it lies at 10 meV). This is the ex- 
pected phonon softening that leads to the high Tc of the 
AlSpGompounds that is evident also in second-derivative 
dataEJ on NbaSn, where the low- Tc material has a large 
peak at 10 meV which disperse to two peaks (one at 6 
and one at 8 meV) in the high- Tc material. Another 
strong indication of enhanced electron-phonon coupling 
is the pronounced softening of the 10 meV shoulder ob- 
served in the phonon density of states of NbaSn on cool- 
ing down from room temperature to 4.2 K as described 
by Schweiss et al.E3 Hence the extracted a^F displays 
the expected phonon-mode softening, with the expected 
energy scales. Another test is to compare the extracted 
a^F to the-|phonon DOS, T(il), measured with neutron 
scatteringLJ. In actuality, neutron scattering weights the 
phonons by the neutron scattering cross section for each 
nucleus forming the generalized DOS G{il.). Nb has a 
cross section approximately twice as big as Sn, but as 
the atomic masses of Nb and Sn are similar, we expect 
both atoms to be in motion for most phonon modes. 
Hence these modes can be excited by interactions of 
neutrons with either Nb or Sn nuclei. This averaging 
effect implies that the phonon density of states F{il) 
and thfi-gejieralized density of states G{il) should agree 
closelyEro for NbaSn. If we take the ratio of the tun- 
neling data to the neutron data, we find a large peak in 
a^(fi) for energies ±>€low 10 meV. This agrees with theo- 
retical calculationiB, which predict a large (fl) for the 
low-energy phonon modes (below 10 meV). 
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Fig. 3, Freericks et al., Phys. Rev. B 

FIG. 3. Eliashberg functions for a tunneling inversion in 
NbaSn for nonconstant density of states (thick) and constant 
density of states (thin). In the nonconstant density of states 
analysis, the spectral weight is enhanced at high energies. 

Next we consider the reduced density of states for both 
cases in Eig. |4| The thick dashed curve is the experimen- 
tal data,El the thin dashed line is the horizontal axis, the 
thick solid line is the nonconstant DOS case and the thin 
solid line is the constant DOS results. There are three im- 
portant points to note about these curves; (i) the fits are 
poor at the lowest energies, (ii) the nonconstant density 
of states produces a more rapid "overswing" and return- 
to-zero at about 40 meV, and (iii) the amplitude of the 
overswing is enhanced in the region from 40 — 60 meV. 
As neither of the fits is particularly good, we cannot con- 
clude from this work that the inclusion of nonconstant 
density of states effects alone is sufficient to explain the 
A15 tunneling data. But we do see that the incorpo- 
ration of a nonconstant density of states definitely pro- 
vides improvements. It is not clear what else is needed. 
Part of the problem may be with the numerical instabil- 
ities of the tunneling-inversion algorithm. Alternatively, 
there may be an intrinsic thin proximity layer that al- 
ways needs to be taken into account regardless of the 
quality of the junction. Other factors that may be im- 
portant include anharmonicity, anisotropy, and impurity 
scattering. It appears unlikely that disorder is the expla- 
nation, since disorder tends to reduce the magnitude of 
the RDOS, not increase it, as is needed. The proximity- 
effect explanation is also hard to support, because the 
constant DOS analysis found the proximity layer to be 
vanishingly thin. If that conclusion holds true for the 
nonconstant DOS analysis as well, then this would not be 
a viable explanation either. The effects of anharmonicity 
at low temperature should be explainable within a quasi- 
harmonic approximation, unless there is a preformed-pair 
phase. Hence, we believe the most likely cause of the 
discrepancy is from anisotropic effects. It is conceivable 
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that tunneling barriers grow differently on different in- 
dexed surfaces. The strong dependence of the tunneling 
conductance on barrier thickness thus could lead to a di- 
rectional selectivity of the tunneling current and thereby 
anisotropy effects would influence the tunneling current 
even in the case of polycrystalline films. Incorporation of 
such anisotropic effects is beyond the scope of this work. 
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FIG. 4. Reduced density of states [from Eq. (|2^)] for non- 
constant (thick) and constant (thin) density of states. The 
thick dashed line is the experimental data (scanned and dig- 
itized) and the thin dashed line is the horizontal axis. Note 
the poor quality of the fit at the lowest energies, and the 
sharper overswing at high energies for the nonconstant den- 
sity of states case. The poor fit arises from the fact that A 
must be increased in order to produce the overswing, but this 
reduces the agreement at low frequencies. The results shown 
represent the best compromise for fitting the entire experi- 
mental spectrum. 
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FIG. 5. Density of states for NbsSn. The solid curve is 
the band-structure density of states as calculated with den- 
sity functional theory. The dashed curve is the quasiparticle 
density of states calculated with the fit a^F. The chemical 
potential lies at tj = 0. Inset shows the full band-structure 
density of states, with the dashed box indicating the region 
blown up in the main figure. The units of the vertical axis for 
both the main plot and the inset are states per spin per meV 
per unit cell. 

In Fig. ^, we plot the real and imaginary parts of x for 
(a) the normal state and (b) the superconducting state. 
This function vanishes for the case of a constant density 
of states. The chemical potentials are = 15.70 meV 
and fis — 15.97 meV in the normal and superconduct- 
ing states. The value of the real part of x is of this 
order of magnitude. The normal and superconducting 
self-energies differ only at the lowest energies. 



We conclude with a discussion of the properties of the 
solutions to the many-body problem. The electronic den- 
sity of states is shown in Fig. |^. The solid curve is 
the band-structure density of states as calculated within 
density-functional theory. The dashed curve is the quasi- 
particle density of states in the normal state at T = 
calculated with the fit a'^F. It depicts the behavior ex- 
pected. The density of states is unchanged at the chem- 
ical potential (w = 0) because the self-energy is mo- 
mentum independent, the peak is narrowed by a factor 
of about 3 (due to "1-l-A" narrowing), and the overall 
density of states is smoothed out due to lifetime effects. 
In the inset, we show the full band-structure density of 
states used in the calculations. The dashed box indicates 
the region shown in the main figure. 
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FIG. 6. Self-energy function x{^) plotted in the (a) nor- 
mal and (b) superconducting states. These functions are 
nearly identical, except for small differences at the lowest en- 
ergies. 

The renormalization function Z{uj) is plotted in Fig. 
for (a) the normal phase and (b) the superconducting 
phase. The solid lines are for the nonconstant density 
of states case and the dashed lines are for the constant 
density of states case. Note how in the nonconstant den- 
sity of states case, the effective strength of the electron 
phonon coupling, measured by Zjv(O) — 1 is closer to 2.1 
than the value of A, which is 2.7. It is the former value 
that is the true measure of the electron-phonon coupling 
strength with a nonconstant density of states. Note that 
the main differences between the nonconstant and con- 
stant density of states calculations is that the overall scale 



is larger for the latter. These functions vary from the 
normal to the superconducting state only at low energies 
as expected. Note further, that the real part of Z can 
dip below 1 for the nonconstant density of states case (it 
never does for a constant density of states). 
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FIG. 7. The renormalization function for (a) the normal 
state and (b) the superconducting state. The solid curves are 
for a nonconstant density of states and the dashed curves are 
for a constant density of states. Note that HeZ can dip below 
1 for the nonconstant density of states case. 

Finally, the superconducting gap function A(cl') is 
shown in Fig. ^. The results for the nonconstant density 
of states (solid line) and the constant density of states 
(dashed line) are nearly identical at both low and high 
energies. In the range from 30 meV to 130 meV, the 
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curves deviate from each other: in the nonconstant den- 
sity of states ease, the peak (real part) and dip (imag- 
inary part) form faster than in the constant density of 
states case. This is what produces the sharper overswing 
in the RDOS for the nonconstant density of states calcu- 
lation. 
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FIG. 8. Superconducting gap function for (solid) noncon- 
stant density of states and (dashed) constant density of states. 
The two curves are nearly identical at low and high energies, 
but differ from about 30 meV to 130 meV. The sharp over- 
swing in the experimental RDOS data is better reproduced 
by the more rapid formation of the high-energy peak (dip) in 
the real (imaginary) part of the nonconstant-density-of-states 
gap function. 



dependence of the density of states is enough to com- 
pletely understand the tunneling data. It may be that 
the tunneling is modified by anisotropic effects, by a nar- 
row proximity-coupled layer or by other effects such as 
anharmonicity. 
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TABLE I. Calculated properties of the Eliashberg function 

extracted from the tunneling inversion. 





inversion 
DOS 


A 


M* 


Win 

[meV] 


A 
[meV] 


[K] 


error 
max. 


error 
r.m.s. 


nonconstant 
constant 


2.738 
2.501 


0.286 
0.210 


7.072 
6.415 


13.652 
11.130 


19 
23 


0.033 
0.034 


0.0048 
0.0055 



IV. CONCLUSIONS 



We have performed a modified McMillan-Rowell tun- 
neling inversion including the effects of the noncon- 
stant electronic density of states near the Fermi level for 
NbaSn. The Eliashberg function that is found by this 
inversion has a larger value of A than generally accepted 
for this material, but the effective value of A derived from 
the renormalization function is more reasonable. Our ex- 
tracted a^F(f2) also has a number of expected features 
including large peaks at about 6 and 8 meV representing 
the soft phonon modes, and a sharply peaked a^{fl) be- 
low 10 mcV which agrees with the theoretical predictions. 
Even taking into account the energy dependence of the 
density of states, we are unable to produce excellent fits 
of the tunneling data, though we can better reproduce 
the overswing observed at high frequencies. We believe 
the problem is partially numerical, as the tunneling in- 
version tries to force sharp spikes into a'^F, but the fit 
is still too poor at the lowest energies and in the over- 
swing region to say that properly including the energy 
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